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Abstract 


The estimation method of the relative density of Culex tritaeniorhynchus summorosus 
larvae in the paddy-fields of an area as a whole was investigated, based on the 
variances of transformed variate y = log (x+1) within and between paddy-fields, 
where x is the number of larvae per dip. As a result, it was shown that only one 
dip should be taken in each paddy-field, and instead, the number of paddy-fields 
examined should be increased, in order to estimate the relative density of larvae in 
the area in maximum efficiency. If one dip is taken in a paddy-field, the estima- 
tion precision becomes higher, of course, when the number of paddy-fields examined 
is increased. In the case of the same number of paddy-fields being examined, the 
higher the mean larval density is, the higher the estimation precision is. 


Introduction 


In a previous paper (Wada et al., 1971), 
it was reported that a sequential sampling 
method is useful to classify the relative 
density of larvae of Culex tritaeniorhyn- 
chus summorosus which is the main vector 
of Japanese encephalitis in Japan, in each 
paddy-field into one of the three density 
levels, low, moderate and high, at a 


given statistical reliability by a relatively 
small number of dips. The sequential 
sampling method can probably be used 
for the classification of the larval density 
also in the whole paddy-fields of an area 
based on the variation in the number of 
larvae between paddy-fields. However, 


not only the classification of the relative 
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density, but also the estimation of the 
relative density itself of summorvosus lar- 
vae in the paddy-fields of an area asa 
whole is often required in ecological stu- 
dies of this mosquito especially in relation 
to the epidemiology of Japanese encepha- 
litis. The estimation method for this 


Materials and 


Materials used to inquire into the me- 
thod of estimating the relative density of 
C. t. summorosus larvae in the whole 
paddy-fields in an area were the same as 
used in the previous paper (Wada et al., 
1971). Mosquito larvae (including pupae) 
were collected by a dipper of 15 cm in 
diameter and 3 cm in depth with a wocden 
handle of 1.2m in length, in paddy 
fields of Mogi near Nagasaki City from 
June to October, 1968. At the outset a 
collector stood at a point on the side of 
paddy-field, and took a dip at the site 
which was thought to be most preferable 


for larvae within the reach of the dipper, 


Then, the collector moved by about 5 m 
along the side and took a dip in the same 
way. In each paddy-field, the dipping 
was made ten times, and the number of 
summoyosus larvae was recorded in each 
dip. The number of paddy-fields thus 
examined was 221 in total. 

The method for estimating the relative 
density of C. t. summorosus larvae or 


the mean number of larvae per dip in 


purpose was inquired into in the present 
paper. 

We are very grateful to Dr. N. Omori, 
the former professor of the Department 
of Medical Zoology, Nagasaki University 
School of Medicine, under whose guidance 


this study was carried out. 


method 


paddy-fields of a particular area in maxi- 
mum efficiency will be discussed, based 
on the variances of the number of larvae 
within and between paddy-fields. As the 
distribution of the number of larvae per 
dip was proved to follow the negative 
binomial distribution with a common +% 
(Wada et al., 1971), the transformation 
of y = log (x + 1) was adopted (Bartlett, 
1947) to normalize the data in the 221 
paddy-fields in each of which ten dips 
were taken, where x is the number of 
larvae in each dip. Firstly the variances 
of y within and between paddy-fields, 
cj and sg, was estimated, and second- 
ly the number of dips within each 
paddy-field that should be taken for 
maximum efficiency with a given cost 
was determined. When the number of 
dips in each paddy-field was thus fixed, 
the relation between the number of 
paddy-fields examined and the estimation 
error of the relative density of larvae in 
the paddy-fields in an area as a whole 


was given. 


Results obtained 


Letting the number of paddy-fields be 


b, the number of dips in a paddy-field a, 


the transformed number of larvae in a 


dip y, the total of y in a paddy-field ż, 
the grand total of y in b paddy-fields T, 
then the correction term, C, and the 
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total sum of squares, SSAB, are calcu- 
lated by 


ab (1) 


SSAB = 3y? — C, (2) 


and the between paddy-fields sum of 
squares, SSB, and the within paddy-fields 
sum of squares, SSA, are given by 


SSB = TE — C; (3) 
SSA = SSAB — SSB. (4) 


The between paddy-fields mean square is 
given by SSB/(b — 1), which is an esti- 


mate of v r + as 2, and the within 


paddy-fields mean square is given by SSA 
/b(a—1), which is an estimate of 


2 where s,” is the population vari- 


A’ B 
ance of between paddy-fields samples and 


o 


g2 


A is that of within paddy-fields sam- 


ples. Thus we have s,” as the estimate 


B 


of as and or as that of oy by the 
following: 


oe SSB/(b—1) a [b(a-1) 


2i = SSA/b(a—1). (6) 


From the present data Sa and s a 


calculated, where b = 221 and a= 10, 


as 


were 


Sy = 0.0116, 
Fi = 0.0152. 


The optimum number of dips per paddy- 
field, x a’ is approximately given (South- 


wood, 1966) by 


C, s? 
=, BOA 
"a -/ ed or ae (7) 

C, s? 
where Cp is the cost of moving to ano- 
ther paddy-field and C A is the cost of 
taking a dip within paddy-field. As it is 
thought to be appropriate that the both 


costs, C, and C are the same, by 


B A’ 


substituting the calculated values of Sy 


and s oF we have 


_ 0.0152 
a 0.0116 


= 1. 


n 


This means that only one dip needs to be 
taken in each paddy-field to estimate the 
mean number of larvae in a particular 
area in maximum efficiency. 


The variance of the mean of y, s$» is 


given by 
2 2 
S S 
$= By 4 (8) 
b b a 


where %, is the number of paddy-fields. 


The number of paddy-fields which must be 
taken depends on the degree of precision 
required. The number of paddy-fields, 


Ny » which needs to be taken to estimate 
the mean of y with the variance of 
2 r l 

d , i. e., the standard error of d, can 


be calculated from 


2 2 
S S 
ng n ni . 


Here, », = 1 and os + vi = (0.0116 
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+ 0.0152 = 0.0268. 

Therefore we have 
0.0268 _ 32 

My (0) 

or 

0.0268 


: t 


Thus we can calculate the values of n, 


with various d, as shown in Fig. 1. It 
is apparent that the number of paddy- 


fields to be examined, Nps decreased 


with the increasing value of d. If we 
wish to estimate the mean of y with the 


accuracy that the standard error, SS 


taken as d is about 0.005, then the num- 
ber of paddy-fields to be examined is as 


many as 1,000, but if s> is about 0.02, 
then 70 paddy-fields will be enough. 
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Fig. {. The relation between the values of the 
level of accuracy (d) and the number 
of paddy-fields to be examined (np, Yy 


obtained from the expression m, = 


0.0268/ d? 
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Fig. 2. The relation between y and D for various np» obtaind from the expression 


n, = 0.0268/ D272, 
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Usually d is expressed as d=Dy, where 
y is the mean of y, and D is a decimal 
which shows the rate of Ss to y. So, 


the expression (11) can be written as 
_ 0.0268 
b 2 ~2 2 
D! y 
From the eXpression (12) we have the 


relation between y and D for various n, >» 


as shown in Fig. 2, and the relation 
between y and ny for various D, as 


shown in Fig. 3. Fig. 2 shows that when 
the same number of paddy-fields (ny) 


are examined, the value of D decreases 
with increasing y, and Fig. 3 shows that 
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when the value of D is the same, the 
number of paddy-fields to be examined 
(n, ) decreases with increasing y. In 


other words, when y is small, D is large 
(the level of accuracy is low), and vice 
versa. For example, in order to estimate 
the mean of y (¥) in the accuracy within 
10% of Sy to y, about 10 paddy-fields are 


enough when y is 0.5, but more than 
1,000 paddy-fields are required when y is 
0.05. 

As Fig. 2 and Fig. 3 were based on 
the values of y which had been transfor- 
med from the number of larvae in a dip 
Cx) by y=log(*+1), it may be convenient 


0.2 0.5 1 2 5 
y (log scale ) 


Fig. 3. The relation y and ny for various D, obtained from the expression 


n, = 0.0268/ D272. 
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to convert back the values of y to the and D for various ny and the relation 


original figures of x by using the relation between x and », for various D were 
y=log (+1) and the relation between x 
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Fig. 4. The relation between x and D for various Bp» obtained by 
converting y into x in Fig. 2. 
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given in Fig. 4 and Fig.5, As x increases 
when y increases, the similar relations 
found in Figs. 2 and 3 are seen in Figs.4 
and 5, namely when * (the mean number 
of larvae per dip) is small the level of 
accuracy is low. Figs. 4 and 5 may be 
used to know the approximate number of 
paddy-fields which need to be examined 
in a given accuracy or to know the level 
of accuracy when a given number of 
paddy-fields are examined, in estimating 
the mean number of larvae per dip in a 
particular area. How many paddy-fields 
must be chosen in planning the larval 
survey of C. ¢. summorosus in paddy- 
fields in a particular area depends on the 
value of D which is determined by the 
purpose of the survey and the abundance 
of larvae. If we take 10% as D, as 
usual, then the approximate number of 
paddy-fields to be examined is 400 when 
the mean number of larvae per dip is 
around 0.2 on an average, 90 when the 
mean is around 0.5 and 30 when the 
mean is around 1.0. In usual breeding 
situation of C. t. summorosus larvae in 
Nagasaki area, 100 paddy-fields may be 
sufficient for most purposes in relation to 
Japanese encephalitis epidemiology. 

As mentioned above, the relative den- 
sity of larvae in water-filled paddy-fields 
can be estimated with a given precision, 
however the total number of C. f. 
summoyosus larvae breeding in a particu- 
lar area as a whole is often required for 
the seasonal changes. The total number 
depends on the density of larvae in the 
water-lodged paddy-fields which are the 
main breeding place for this mosquito 
larvae, and also on the surface area of 
the whole paddy-fields in the area. 


Accordingly, for the index of the total 
number of larvae, it is advisable to 
multiply the mean number of larvae per 
dip by the rate of water-lodged surface 
area in the paddy-fields in the area. 

In the present paper, the number of 
larvae per dip (+) was transformed into 
y= log (+1) so that the variance of the 
transformed variate was independent of 
the level of mean, and then the estima- 
tion method of relative density of the 
larvae was discussed. Here, the relation 
between the variances before and after 
transformation will be mentioned. 

If a variate x is transformed into y by 
y=f(x), then 

V=fOM) +f (maem) bree > 43 
where m is the expectation of x. There- 
fore, approximately 

(y=M)=f (m)(x—-m), E) 
where M is the expectation of y. Letting 


the variance of y be 7” and that of x be 


2 


a”, we have 


o = {pm}? 0? ) 


If a is related to m as 2 = 7 m (A 


2 becomes nearly cons- 


is a constant), 7 
tant by adopting the transformation y= 
log(++1), and in negative binomial distri- 
bution this transformation seems likely to 
prove good enough in many cases (Beall, 
1942). Since f (4) =loge / (+1), f (m) 
=0.4343 /(m+1), and accordingly from 
expression (5), 
a” = 0,189 a? / (m4+1)" . (6) 

Thus, we have the approximate relation 


between the variances before and after 
the transformation, 7? and o? . If m is 


fairly large, then 2’? comes near to 0.189 


: 2 2 
» since =) m. 
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